Relativistic ionization-rescattering with tailored laser pulses 
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The interaction of relativistically strong tailored laser pulses with an atomic system is considered. 
Due to a special tailoring of the laser pulse, the suppression of the relativistic drift of the ionized 
electron and a dramatic enhancement of the rescattering probability is shown to be achievable. The 
high harmonic generation rate in the relativistic regime is calculated and shown to be increased by 
several orders of magnitude compared to the case of conventional laser pulses. The energies of the 
revisiting electron at the atomic core can approach the MeV domain, thus rendering hard x-ray 
harmonics and nuclear reactions with single atoms feasible. 
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In intense laser-atom interaction phenomena the 
rescattering concept plays a central role 0. Here the 
electron is ionized, propagated in the continuum by 
the laser field, and finally driven back and scattered at 
the ionic core in the case of above-threshold ionization 
(ATI) 0|, ionizing further bound electrons in the case of 
non-sequential double ionization (NSDI) 3], or recom- 
bining with the ionic core with high-harmonic genera- 
tion (HHG) Q. It is highly desirable to increase the 
rescattering electron energy [5l as it can be employed 
to generate higher harmonics la, to image attosecond 
dynamics of nuclear processes p], or to initiate nuclear 
reactions with a single atom/molecule The elec- 

tron energy increase can not be achieved by a straight- 
forward increase of the laser intensity. When the laser in- 
tensity approaches the relativistic regime, the laser mag- 
netic field effect starts to play a role by inducing a drift 
of the ionized electron in the laser propagation direction 
which severely suppresses the probability of the electron 
revisiting the ionic core. There are several attempts to 
circumvent this effect in order to increase the efficiency 
of rescattering, particularly, using relativistic ions which 
propagate in laser propagation direction 0, Q, or us- 
ing two counter-propagating laser beams with linear 
or equally handed circular polarization , or generating 
harmonics via exotic positronium atoms in strong laser 
fields 0]. O n a different front, new techniques for gener- 
ating attosecond puls es have recently emerged based on 
HHG in gases [r3.ll3| and in plasmas interaction 
The usage of attosecond pulse trains (APT's) to enhance 
HHG in a strong laser field has been shown in 

In this letter we show how efficient recollisions in the 
relativistic regime are feasible by employing specially tai- 
lored strong laser pulses in the form of APT's (see FigUJ. 
As an indicator of strong field recollision phenomena we 
consider the HHG process from an ion in a relativistically 
strong tailored laser pulse (see FigEJ). While in 0] the 
weak APT serves to solely control the initial conditions of 
the ionized electron in the strong laser field, we consider 
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FIG. 1: By using optimized tailored pulses in the form of an 
APT instead of conventional sinusoidal pulses it is possible to 
reduce significantly the drift of the ionized electrons due to 
the Lorentz force of the laser field (the ionized electron tra- 
jectory in the sinusoidal pulse is shown as a dashed line, in 
the tailored pulse as a solid line). Therefore, the rescattering 
probability can be increased by several orders of magnitude 
in the relativistic regime of laser-atom interactions, x and z 
are the laser polarization and propagation directions, respec- 
tively. 



HHG in relativistically strong APT. The temporal tailor- 
ing of the laser pulse is aimed to concentrate the ionizing 
and afterwards the accelerating laser forces in short time 
intervals within the laser period (maintaining the aver- 
age intensity of the pulse constant) which substantially 
reduces the relativistic drift. This results from the fact 
that in the tailored laser pulse, fragments in the electron 
trajectory are avoided, in contrast to the sinusoidal laser 
pulse, in which the electron acceleration is compensated 
by deceleration without a net energy gain by the elec- 
tron, but during which nevertheless the electron would 
continue to drift in the laser propagation direction (see 
Fig. EJii))- This substantially decreases the time span 
during which the electron moves with a relativistic ve- 
locity. The latter results in a shorter drift in the laser 
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FIG. 2: The tailored laser pulse in the form of an APT, with 
each pulse duration in the train being r = 0.127T a.u.: (i) 
Phase dependence of the electric field (Eirj)) of the tailored 
pulse (solid line) compared with a sinusoidal field (dashed 
line) with equal average intensities and central angular fre- 
quency w — 0.1 a.u.; (ii) Frequency spectrum (S(k)) of the 
tailored pulse. The frequency components are phase locked. 



propagation direction, and therefore, in an increase of 
the recombination probability in the case of the tailored 
pulse. 

The scenario of the electron motion in the tailored laser 
pulse looks as follows. The trajectory of the rescattering 
electron begins at the end of the period with maximal 
field strength of the laser pulse, shortly before the field 
free region. The initial spike of the laser pulse ionizes 
the electron and initiates its excursion with a very low 
initial velocity. The acceleration of the electron takes 
place during the next spike after the field free region. 
One can see from Fig. |3] that in the tailored wave the 
electron gains most of its energy only shortly before the 
recombination and the trajectory of the electron in the 
tailored pulse is spatially significantly shorter than in the 
sinusoidal wave. 

The procedure for tailoring the laser pulse is the follow- 
ing. The laser pulse with temporal form as shown in Fig|3 
has been chosen with the additional condition of vanish- 
ing the field average of the pulse < E > = 0. Then, the 
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FIG. 3: Trajectory and kinetic energy of an HHG electron in 
the tailored wave (solid line) and in a sinusoidal wave (dashed 
line) with an averaged laser intensity of 5 x 10 19 W/cm 2 . (a; 
and z are the polarization and propagation direction of the 
laser field, respectively.) The electron is ionized at a phase of 
-1.24 in the sinus wave and recombines at a phase of 2.68. In 
the tailored wave the electron is ionized at a phase of -2.78 
and recombines at a phase of 3.14. 



Fourier-expansion of this pulse is used via K = 100 har- 
monic terms to optimize the field expansion coefficients 
by maximizing the rescattering (HHG) probability. The 
optimized tailored pulse has the form of an APT which 
includes three key features. The first is a large ionization 
peak with a strong decay after the point of ionization 
(the stronger decay the larger HHG efficiency) . This fol- 
lows by a long region with an only weak field strength 
(the longer the better) which is resumed by a short and 
strong pulse that induces most of the energy of the ion- 
ized electron. The optimized pulse along with its spec- 
trum is shown in Fig. [5] The rescattering rate resulting 
from this pulse can be enhanced by using more frequency 
components in the APT. In the calculated spectra that 
will follow, we employ 100 frequency components. When 
using only K — 30 harmonics in the tailored pulse, the 
HHG yield decreases by approximately 2 orders of magni- 
tude at laser intensity of 5 x 10 18 W/cm 2 and by 5 orders 
of magnitude for 5 x 10 19 W/cm 2 . A possible experimen- 
tal way for the above mentioned tailoring of the laser 
pulses is to employ relativistic harmonics generated dur- 
ing a laser-foil interaction in the so-called sliding mirror 
regime |l5j which yields very strong APT's with a con- 
version efficiency of a few percent. Additionally the tech- 
nique for the spectral filtering should be applied to 
form the necessary spectral dips as shown in Fig. [21 Sig- 
nificantly less spectral filtering would be necessary when 
using K — 30 harmonics and a APT spectrum with one 
hump, though in this case the HHG yield will be reduced 
as indicated above. 

Given the tailored laser pulse, we proceed to calculate 
the HHG intensity. Our investigation of HHG in an opti- 
mized tailored laser pulse is based on the solution of the 
Klein-Gordon equation in the strong field approximation 
(SFA) |l7|. We consider the HHG process of an atomic 
system by a linearly polarized plane laser field for the 
relativistic parameter regime. The direction of the laser 
field polarization is along the x axis, that of the magnetic 
field along the y, and of the laser propagation along the 
z axis. The amplitude of HHG within the SFA based on 
the Klein-Gordon equation in the single-active electron 
approximation is given by the following expression |Tsj 
(atomic units are used throughout the paper): 

M n = -if dV J d 4 x"{$(x'y 

x V H (x')G v {x',x")V A {x")<$>{x")} (1) 

with n the harmonic order, $ the bound state 
wave function, Vh(x) = 2Ah{x) ■ (p + A(r])/c), 
Va(x) = 2iV(x)/c 2 d t + V{x) 2 /c 2 , V(x) the atomic po- 
tential, c the speed of light, p the momentum operator of 
the electron, A (77) the vector potential of the laser field 
in the radiation gauge, k^ — (u)/c, 0, 0, k) its 4-wave- 
vector, 77 = k^x^, x^ = (ct, x) the time-space coordinate, 
A-h(x) the matrix element of the vector potential of 



the high harmonic field in the second quantization for 
an one photon emission process, and G v (x',x") the 
Klein-Gordon Volkov Green function given in Ref. 
$ is an eigenstate of the physical energy operator in the 

radiation gauge and can be approximated as <fr(x) = 

0o(x) exp {-i[{c 2 - I p )t + x • A/c]} /y/2(c? ' - I p ) (H, 
with the nonrelativistic ground state wave function 
4>q (x). Introducing the integration variable transforma- 
tions Tj' = k^x' and 77" = k^x"^ and taking into account 
that the atomic wave function characteristic length is 
much smaller than the laser wavelength, the amplitude 
of HHG can be written: 

/OO rr)' r 

drf dr," / d 3 qm ff (q,,',,") 
-OO J — OO J 

xex P {-i[S(q,t ? ',7/ / )-nT/]} (2) 

with the HHG matrix element for emission in laser po- 
larization direction 



h 1 , f/\ c 2 { Px +A{j)/c) 
m {p,V >V ) = 9 

A(tj') h 



(3) 



x 0* p 



A(V) 



C 2 )|V|0 



with the Fourier transform of the ground state wave 
function </> and e p = ^/c 4 + c 2 p 2 . £(p, 77, ?/) — 
f^,df](e p (ri) — c 2 +Ip)/Lu is the quasi-classical action 
and the relativistic energy is given by 



/c • p 



(p + A(r?)/2c) • A{t))/c. 



The differential rate of HHG is then given by 
dw n 



dn 



MJ 



(4) 



(5) 



with the solid angle of emission fl. 

Using the low-frequency approximation Klu <C I p with 
the ionization potential I p , the integral in Eq.Q can be 
calculated via the saddle-point method. The saddle point 
conditions along with the cutoff condition determine the 
ionization and recombination phases rj" , fj', as well as the 
cutoff harmonic number n: 



£q(?/) = c 2 -Ip 
Hv") = c 2 -I p 



drf 



= 0, 



(6) 



with the drift momentum of the ionized electron q x — 

J?!,dfjA(fj)/c . n _ 



rj'-i)" 



% = 0, q z 



= , with q 2 m = 



-(1/c 2 ) p'„ dfjA{f)) 2 /(f)' -i)") 0. In tailoring the laser 
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FIG. 4: Harmonic emission rate in laser polarization direction 
via log 10 (dw„/dQ,) in Eq. 0, as function of the harmonic 
energy with a laser intensity of (i) 5 x 10 18 W/cm 2 and (ii) 
5 x 10 19 W/cm 2 . The main angular frequency is uj — 0.1 a.u.. 
The ionization potential is (i) I p = 28 a.u. and (ii) I p — 62 
a.u., respectively: (a, grey) within the dipole approximation 
and a sinusoidal field, (b, dashed) with respect to the Klein- 
Gordon equation and a sinusoidal field, (c, black) with respect 
to the Klein-Gordon equation and a tailored pulse. 



pulse, we fix the harmonic cutoff h and by changing the 
shape of the pulse, minimize the imaginary part of the 
ionization time Im{fi"} which determines the plateau 
height and, therefore, the rescattering efficiency. Fur- 
ther, the atomic potential V(x) is approximated by a 
zero-range potential |20| . The ionization potential I p is 
adapted to the threshold of the barrier suppression ion- 
ization model. 

The HHG spectrum with the tailored pulse is com- 
pared with two calculations for the harmonic emission 
using a sinusoidal field. While the first calculation is 
based on the dipole approximation (DA), the second one 
is fully-relativistic. We consider two regimes: the moder- 
ately relativistic regime with a laser intensity of the tai- 
lored pulse 5 x 10 18 W/cm 2 and the strongly relativistic 
regime with a laser intensity of 5 x 10 19 W/cm 2 , each with 
a laser ground angular frequency of 0.1 a.u., correspond- 
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ing to a laser wavelength of A = 456 nm. Taking into 
account a few percent efficiency in the APT generation 
via laser-foil interaction 0] , as well as about 20% inten- 
sity reduction because of spectral filtering in the case of 
K = fOO spectral components in APT, the initial laser 
intensity before tailoring should be three orders of mag- 
nitude larger. For pulses involving fewer high harmonics 
this loss may be reduced. 

In the moderately relativistic regime (see FigQJi)), the 
harmonic emission rate with conventional laser field is 
suppressed by over seven orders of magnitude in the 
relativistic treatment compared with the results in the 
DA. Whereas the tailored pulse yields a harmonic emis- 
sion rate that is nearly four orders of magnitude higher 
than the rate in the dipole approximation. Further the 
spectrum induced by the tailored pulse as well as the 
spectrum induced by the sinusoidal field using the DA is 
highly oscillating, resulting from interference of at least 
two trajectories of the ionized electron. In the case of 
the tailored pulse, there are exactly two, since the tra- 
jectories with multiple return to the ionic core are very 
unlikely. 

The HHG spectrum in the strong relativistic regime 
with a laser intensity of 5 x f0 19 W/cm 2 is shown in 
Fig. IHii). Whereas the tailored pulse spectrum is still 
more intense than the spectrum in the DA with a con- 
ventional laser field, the relativistic treatment of HHG 
with sinusoidal fields yields a rate that is negligibly small. 
From the HHG spectrum one can see that cutoff energies 
of more than 1 MeV are reachable with the probability 
not less than in the DA. 

The enhancement of the HHG rate via tailored pulses 
compared with the one via a sinusoidal field is due to 
the suppressed relativistic drift. This can be deduced by 
comparing the trajectory of the ionized cutoff electron 
driven by the tailored pulse with one driven by the sinu- 
soidal pulse, see Fig. 034). The excursion of the ionized 
electron in polarization and propagation direction in the 
tailored pulse is shorter, though the averaged intensity 
of the pulses are the same. Moreover, the HHG emission 
angle with respect to the laser polarization direction in 
the tailored pulse is significantly smaller than in the con- 
ventional sinusoidal field (see Fig. Et 1 ))- ln other terms, 
the suppressed drift can be deduced by the value of ini- 
tial electron velocity in the laser propagation direction, 
necessary for the electron rescattering in the relativistic 
regime. In the tailored pulse it is much smaller compared 
to that for a sinusoidal pulse. In fact, in the strongly 
relativistic regime, the numerical values for the initial 
velocity in the tailored pulse is -6 a.u., whereas for the 
sinusoidal field -97 a.u.. The initial velocity of the tun- 
neled electron determines the exponential damping fac- 
tor of the electron tunneling probability, yielding a higher 
ionization rate and, therefore, a higher HHG efficiency, 
in the tailored pulse as compared to the sinusoidal field. 

The reason for the increased HHG rate with the tai- 



lored laser pulse over the dipole approximation rate is 
that the laser electric field strength at the time of ioniza- 
tion is higher in the tailored pulse at the same average 
laser intensity. The numerical values are \E(f}")\ = 74 
a.u. for the tailored pulse and \E(fj")\ = 36 a.u. for the 
sinusoidal field for a laser intensity of 5 x 10 19 W/cm 2 . 

Further, one can make a statement about the stability 
of the HHG process driven by such a tailored laser pulse. 
A random variation of the amplitude of the frequency 
components by 5 percent, and of the phases within the 
range not exceeding 7r/100, results in a variation of the 
HHG yield within one order of magnitude. 

Concluding, employing specially tailored laser pulses 
in the form of an APT allows to suppress the relativis- 
tic drift of the ionized electron and to realize ionization- 
rescattering in the relativistic regime. HHG in the hard 
x-ray domain and the initiation of nuclear reactions with 
single ions thus become feasible. 
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